The even-odd effect in short antiferromagnetic Heisenberg chains 
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Motivated by recent experiments on chemically synthesized magnetic molecular chains we inves- 
tigate the lowest lying energy band of short spin-s antiferromagnetic Heisenberg chains focusing on 
effects of open boundaries. By numerical diagonalization we find that the Lande pattern in the 
energy levels, i.e. E(S) oc S(S + 1) for total spin S, known from e.g. ring-shaped nanomagnets, can 
be recovered in odd-membered chains while strong deviations are found for the lowest excitations in 
chains with an even number of sites. This particular even-odd effect in the short Heisenberg chains 
cannot be explained by simple effective Hamiltonians and symmetry arguments. We go beyond 
these approaches, taking into account quantum fluctuations by means of a path integral description 
and the valence bond basis, but the resulting quantum edge-spin picture which is known to work 
well for long chains does not agree with the numerical results for short chains and cannot explain 
the even-odd effect. Instead, by analyzing also the classical chain model, we show that spatial fluc- 
tuations dominate the physical behavior in short chains, with length N < e vs , for any spin s. Such 
short chains are found to display a unique behavior, which is not related to the thermodynamic 
limit and cannot be described well by theories developed for this regime. 

PACS numbers: 75.50.Xx, 75.10.Jm 



INTRODUCTION 



The antiferromagnetic Heisenberg model is the appro- 
priate starting point to understand magnetism in a large 
variety of different materials where strong electron cor- 
relations are important. Within its context the quantum 
nature of magnetism reveals itself, leading to a vast ar- 
ray of fascinating phenomena. Accordingly this model 
has been the topic of numerous experimental and theo- 
retical works in different combinations of magnetic lat- 
tices and spin magnitudes. Indeed, many interesting 
theoretical concepts have been developed specifically for 
the antiferromagnetic Heisenberg model in the thermody- 
namic limit r — including spin waves, effective Hamilto- 
nians, quantum field theories, and hydrodynamic meth- 
ods. However, as will be shown in this paper most of 
these methods fail when describing finite (but possibly 
large) clusters of spins, which show a unique behavior 
that strongly depends on the boundary conditions and 
topology and therefore form a separate class of magnetic 
materials. 

The properties of the antiferromagnetic Heisenberg 
model on small clusters such as dimers, trimers and 
tetramers, which are relatively simple to understand, 
were mainly of relevance to chemists, as the Heisenberg 
exchange describes magnetic interactions between metal 
ions in polynuclear metal complexes very well,— and also 
allowed the investigation of the origin of the fundamen- 
tal magnetic interactions^ However, the advances in the 



synthesis of metal complexes during the last fifteen years 
have produced nanosized magnetic molecules with up to 
a few dozens of magnetic metal ions interacting with each 
other via Heisenberg exchange, creating the new class of 
molecular nanomagnets,- 1 ^ Molecular nanomagnets be- 
long to the mesoscopic regime, and their magnetism can 
be considerably more complex than in few-membered 
metal clusters and possess a quantum many-body char- 
acter, yet they may not contain enough spin centers to be 
appropriately described by the methods and techniques 
developed for extended systems. Apart from synthetic 
chemistry, artificial engineering of quantum spin clus- 
ters has also emerged in recent years J^~— Here clusters of 
magnetic ions have been fabricated directly on insulating 
surfaces, and their magnetic properties were measured 
with scanning tunneling microscopy. 

Nanosized spin clusters are ideal to study basic ques- 
tions of quantum mechanics in mesoscopic systems, such 
as the efficiency of different metal centers or topologies 
towards a desired magnetic property, or the transition 
from the quantum to the classical regime for larger ion 
spins They also provide an ideal testing ground for 
the validity of different theoretical models. In this paper 
we consider the antiferromagnetic Heisenberg model with 
a small number of spins of magnitude s placed along a 
chain of length N with open ends, 
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where denotes a spin-s operator on site i and J > for 
antiferromagnetic interactions. We will show that despite 
the apparent simplicity of the model the detailed aspects 
of the spectrum are highly non-trivial. 

The model for an infinite chain or one-dimensional an- 
tiferromagnetic Heisenberg chain (AFHC) has attracted 
enormous interest, especially after Haldane's conjecture 
of a fundamental difference between integer and half- 
integer spin chains, which is by now well established^ 2 - - — 
Finite but long chains have also received significant at- 
tention, but mainly for connecting numerical results to 
the thermodynamic limit via scaling^ or for understand- 
ing boundary^ impurity^ and doping effects. 1 — On the 
other hand, very little has been done for short chains 
with s > 1/2, in which the Haldane gap is always smaller 
than the finite-size excitation energy, and the difference 
between integer and half-integer spins is thus expected 
to become irrelevant. 

Short AFHCs have recently become accessible exper- 
imentally, e.g., through the synthesis of molecular Cr 3+ 
(s = 3/2) modified wheels and horseshoes J^r— Two Fe 3+ 
containing ring molecules were also studied, which mag- 
netically represent s = 5/2 chains of lengths N = 7 
and 17i 26 i 27 Magnetic and inelastic neutron scattering 
measurements on the Crg and Crr horseshoes (s = 3/2, 
N = 6 and 7) have provided a detailed view on the 
magnetic excitation spectrum of these clusters ; 1 ^ 20 and 
pointed to systematic differences between chains with 
even and odd number of spins. The expected even-odd 
effects were found, such as total spin 5 = and 5 = s 
in the ground state for even and odd N respectively, 
which can obviously be associated to the interplay of 
the open boundary conditions and the symmetry require- 
ments. However, a more subtle but striking difference in 
the energy spectrum has also been noted in addition, 
which we here refer to as the even-odd effect. 

This even-odd effect manifests itself in the lowest en- 
ergy band as a function of total spin 5, which depending 
on the context is known as the tower of states, quasi- 
degenerate joint states, or the L bandit— In this band 
the energies are expected to increase as E(S) oc 5(5 + 1), 
and for bipartite spin systems it is generally possible to 
use a simplified model, the Hab Hamiltonian, to de- 
scribe it. The Hab model has been very successfully 
applied to even-membered antiferromagnetic Heisenberg 
rings (odd-membered rings will not be considered in this 
work)^ and experimentally confirmed in fine detail for 
the Cr§, CsFes and Feis molecular wheels^Sr— Further 
examples are the Mn-[3x3] grid and the Fe3o Keplerate 
molecules i^Sr— For short chains with open ends, how- 
ever, the E(S) oc 5(5 + 1) approximation of the L band 
appears to work very well for odd chains, while for even 
chains it surprisingly fails in the low energy sector. 

We here explore alternative approaches to extract the 
essential physics of the AFHC model for short chains. 
In particular, we carefully reanalyze the Hab model, 
and compare numerical diagonalization results for Eq. (TT]) 
with the predictions of field theoretical approaches, the 



classical AFHC, and parent valence-bond Hamiltonians. 
These approaches differ in their treatment of the quan- 
tum and spatial fluctuations and allow us to get insight 
into the roles played by them. 

As main results we will demonstrate that the edge-spin 
picture, which has been firmly established for long chains 
for s = liSri 3 . and higher g ) 15 i 44 i 45 anc j was conjectured to 
stay robust up to very small iV ) 15 i 45 is not consistent with 
our exact numerical data, justifying our distinction into 
short and long chains. Furthermore, for spins s > 3/2 
the even-odd effect has converged already closely to the 
classical behavior, yet quantum fluctuations arc still no- 
ticeable. We finally conclude that short AFHCs posses 
simultaneously a classical and a quantum character, blur- 
ring the limit between classical and quantum behavior. 

The paper is organized as follows: in Sec. |TT]the even- 
odd effect is introduced, and the Hab model is analyzed 
in Sec. Mil Section IIVI deals with the predictions of the 
0(3) nonlinear a- model, which originates in a path in- 
tegral representation. Section [V] goes away from semi- 
classical descriptions to investigate quantum mechani- 
cal effects on the extreme quantum s — 1/2 case. The 
classical limit of the AFHC is investigated in Sec. IVI| 
and spin density and correlation effects in Sec. I VIII The 
valence-bond model description of the AFHC is analyzed 
in Sec. IVIlJ and finally, Sec. IIXI presents the conclusions. 
The Appendices that then follow include relevant infor- 
mation on various topics of the main text. 



II. THE EVEN-ODD EFFECT 

The different symmetry properties between even- and 
odd-membered AFHCs naturally generate differences in 
the magnetic properties, which show up in various quan- 
tities giving rise to various "even-odd" effects. The most 
obvious difference is in the ground state, which has to- 
tal spin S g — for even and S g = s for odd chains due 
to the residual spin, as can be inferred from the clas- 
sical antiferromagnetic spin configuration in the ground 
state, or the theorem of Lieb and Mattis4£ Pronounced 
differences are also present for example in the spin den- 
sity and spin correlation functions (see Sec. IVIII) . which 
are however largely dictated by symmetry and boundary 
considerations. 

The even-odd effect observed in this paper is related 
to pronounced deviations from an approximate E(S) oc 
5(5 + 1) energy dependence, and through this is asso- 
ciated to the Hab model, or more generally the L- and 
_E-band picture of the excitations in not too large bi- 
partite spin clusters (the _B-band collectively denotes the 
next higher lying rotational bands). 7 In this section we 
define and characterize the even-odd effect. The Hab 
model and the L&i?-band picture will be described and 
reanalyzed in more detail in the following section. 

The Hamiltonian of Eq. ((T|) is SU(2) invariant and the 
eigenstates are organized as multiplets of 25 + 1 degen- 
erate states according to 5. We focus on the lowest spin 
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FIG. 1: (Color online) Energy spectra of antiferromagnetic Heisenberg chains as function of S(S + 1) for (a) N = 6, s = 3/2, 
(b) N = 7, s = 3/2, (c) N = 6, s — 3, and (d) N — 7, s = 3. Panels (a) and (b) show the full spectrum, while (c) and (d) the 
low lying spectrum. In panels (a) and (b), the solid lines indicate the slopes §af (red line), gF (green line), and the prediction 
of the Hab model (blue line), which establishes also an upper bound for the minimal states in each S sector. The L band is 
clearly visible in all panels; in panels (c) and (d) it is emphasized together with the E band by the thick underlying bars. 



multiplet in each 5 sector, i.e., the L band. Starting 
from the ground state, the energies of the lowest spin 
multiplets increase approximately as 5(5 + 1) (consis- 
tent with the Lieb-Mattis ordering of energies 4 ^) . They 
thus fall approximately on a straight line in an energy 
vs. 5(5+1) plot, which is the L band. The ground state 
spin will be denoted by S g , and that of the state with 
maximal spin or ferromagnetic state by 5/ = Ns. The 
region in the energy spectrum with small (large) values 
of 5 will be called the antiferromagnetic (ferromagnetic) 
region. Energy values are always presented with respect 
to the ground-state energy. The eigenvalues and eigen- 
states of Eq. |T]) have been calculated by full exact, sub- 
space iteration and/or Lanczos diagonalization, with the 
Hamiltonian matrix represented in the product space or 
by irreducible tensor operators^, and spatial symmetry 
was also use d 47 ' 48 (see Appendix lAl for more details). 

For illustration, Fig. [1] shows spectra of AFHCs with 
N = 6 and 7 for both s = 3/2 and 3. While it is clearly 
possible to identify the L band, there arc systematic de- 
viations. In particular, for N = 6 the L band is no longer 
well approximated by an 5(5 + 1) dependence for small 
5. In contrast, for N = 7 the 5(5 + 1) behavior is sur- 
prisingly well obeyed. Depending on the view point, the 
question hence arises why the L band deviates from the 



5(5 + 1) behavior so strongly in the even chain, or why 
is it so well realized in the odd chain. 

In order to characterize the differences between even 
and odd N we consider the (normalized) slope of the L 
band in the E vs. 5(5 + 1) representation as a function 
of 5, 



g(S) 



2 E(S)-E(S-1) 



*AB 



25 



(2) 



with 5 > S g , which is the discretized version of g(S) = 
-£^dE/d[S{S + 1)]. The slope is normalized with the 
quantity A^s/2, which is given for even and odd chains, 
and rings (for comparison), as 

even chain: Aab 
odd chain: 

i\ -t- i 
4J 

even ring: A AB = — . (3) 



AJ(N-l) 

4J 
N+l' 
AJ 

W 



The normalization allows us to directly compare with the 
predictions of the Hab model, for which g(S) = 1 (see 
Sec. IIII[) . It is noted that the difference of for the 

even ring and the chains is of order 0(l/iV 2 ) while that 
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for the even and odd chain is of order 0(1/N 3 ), and that 
the slope emphasizes differences in the antifcrromagnetic 
region. 

Experimentally, the energies E(S) are directly con- 
nected to the low-temperature magnetization curve 
M(B), and their differences can be extracted from it (B 
is the magnetic field). In a magnetic field the energies 
are modified as Eb{S,B) = E(S) — 2/jbBS (assuming 
a gyromagnetic factor equal to 2), and at zero temper- 
ature M is discontinuous with steps of height 2fis at 
fields B s = [E(S) - E(S - l)]/{2f iB ). An energy depen- 
dence E{S) = ^AS(S + 1) with an appropriate gap A 
yields steps at regularly spaced fields Bs = -^^S, and 
deviations from the S(S + 1) energy dependence or a 
constant slope g(S) are detected as deviations from this 
regular field pattern. If one extrapolates through the 
magnetization steps or measures the magnetization at el- 
evated temperatures, the M(B) curve increases linearly 
with field (except at very low fields) , and deviations from 
the E(S) cx S(S + 1) behavior are observed as a field- 
dependent slope or non-constant susceptibility M/B. For 
the classical case s — > oo the slope can in fact rigorously 
be shown to be proportional to the inverse susceptibility 
(Sec. Eg. The energy differences E(S) - E(S - 1) can 
also be measured directly by inelastic neutron scatter- 
ing due to the selection rule \S — S'\ = 1, permitting a 
spectroscopic determination of the slope g(S). 

The slope g(S) is plotted for N = 6 and 7 for s ranging 
from 1/2 to 7/2 or 3 in Fig. [3] For N — 7, which rep- 
resents odd chains, g(S) varies weakly with S (at least 
for s > 3/2), and is close to 1. In contrast, there is a 
strong reduction of the slope for N = 6, the representa- 
tive for even chains, for S < 2s. In comparison to the 
maximal total spin Sf the region of strong deviation is 
thus given by S/Sf < 2/N, and the disagreement with 
the S(S + 1) spectrum does not alleviate with increasing 
s, contrary to naive expectation. The deviation is pro- 
nounced for a wide range of S values for short chains, but 
for long chains where s <C Sf it becomes very limited in 
range and is less important. These observations are a ma- 
jor result of this paper and establish the even-odd effect 
considered here. 

The even-odd effect is most pronounced in the antifcr- 
romagnetic region and hence also the slopes 



9af = g(S g + l), 
9f = g(Sf) 



(4) 



are considered, which characterize the antiferromagnetic 
and ferromagnetic parts of the L band respectively [see 
Figs da), (b)]. In Fig. [3] is plotted as function of N 
for various s. Since Aab in Eq. © is slightly smaller for 
even than for odd chains, the scaling by Aab mitigates 
the even-odd difference in gAF as compared to the non- 
normalized slope. Nevertheless, the difference between 
even and odd chains is obvious. Interestingly, the differ- 
ence between the even and odd gAF grows with N. This 
divergence suggests that the excitations on the L band 
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FIG. 2: (Color online) Slope g(S) as function of S for S > S g 
for (a) TV = 6 and (b) N — 7 chains for different spins s 
(symbols). The horizontal (red) line is the slope g(S) = 1 of 
the Hab model. The arrow and the dotted line in (a) point 
to the value of the slope for 5* = 2s. Lines are guides to the 
eye. 



in the antiferromagnetic region have different physical 
origin for even and odd chains. For much larger N the 
behavior known for long chains will be approached. The 
ferromagnetic slope gF can be calculated exactly from 
the lowest one-magnon energy as^ 
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Hence, as gp is just the classical one-magnon energy for 
a particular wave vector, the even-odd effect is absent in 
the ferromagnetic regime. We note that using Eq. ([S]) 
one finds gF > 1 for N > 3. This is an important differ- 
ence to even-membered rings, where due to translational 
invariance the eigenvalues and eigenfunctions of the low- 
est one-magnon state arc equivalent to those of Hab and 
hence g r p ng = 1 for all ring sizes^i The absence of the 
even-odd effect in the ferromagnetic region is also seen 
in the spin density in the lowest one-magnon state, given 
by 



\Si) = s- — cos | ;t- 



2i- 1 



27V 



(6) 



There is no qualitative difference between even and odd 
N here. 
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FIG. 3: (Color online) Slope qaf for even (blue lines) and 
odd (green dashed lines) chains as a function of the length 
N for different spins s (symbols). The horizontal (red) line 
represents the slope g(S) = 1 predicted from the Hab model. 
Lines are guides to the eye. 




1/2 1 3/2 2 5/2 3 7/2 4 9/2 5 
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FIG. 4: (Color online) Normalized excitation energies 5e{s) of 
the lowest energy states in the spin sectors S < S g for the N = 
7 chain and spins s ranging from 1 to 5 (symbols). For each 
chain an approximate linear dependence on S is observed. 
Lines are guides to the eye. 



TABLE I: Symmetry properties of the lowest spin multiplets 
in each S sector for the AFHC. 



N 


s 


mirror 


spin-flip 






parity 


parity 


even 


integer or half-integer 




/-QS+iVs 


odd 


integer s, S > S g 


+ 1 


(_!)*-. 


odd 


integer s, S < S g 


(_!)*-. 


(_!)*- 


odd 


half-integer s, S > S g 


+1 




odd 


half-integer s, S < S g 







Even-membered chains on the other hand are mirror sym- 
metric about a link and cannot show any local magne- 
tization in the ground state. In this case an alternating 
quantum dimerization is more suggestive if s is not very 
large, where neighboring spins are more strongly corre- 
lated on odd links and less on even links. Another dif- 
ference is that possible residual spin degrees of freedom 
at the edges, the so-called edge spins, first introduced 
in the context of valence bond ground states and long 
s = 1 AFHCs , 40 i 43 ' 49 i 50 are effectively coupled ferromag- 
netically for odd N and antiferromagnetically for even N 
(Sec. IIVC|) . Notably, rings with periodic boundary con- 
ditions always possess mirror symmetry about a link and 
a site for both even and odd N, so that there is no fun- 
damental symmetry difference as in chains. For chains 
the symmetry properties of the lowest multiplets in each 
sector 5, which for 5 > S g represent the L band, are 
summarized in Table HI — in agreement with the findings 
for the s = 1/2 chain in Ref. [El]]. It is also noted that 
multiplets in a specific 5 sector share the spin-flip parity 
of the corresponding state of the L band, whenever this 
is a good quantum number. 



III. THE H AB MODEL 



It is also of interest to consider the lowest energies 
E(S) in odd chains for S < S g = s, since the Hab model 
(Sec. IIIip and the edge-spin picture (Sec. IIV[) predict 
distinctly different trends with 5 for these states. The 
energies in this region can be represented by the normal- 
ized excitation energies 



6e(S) = 



2 [E(S) - E(S g )} 
A AB (N-l)s 



(7) 



shown in Fig. @]for the N = 7 chain for s = 1 to 5. The 
energies get smaller with increasing S, E(S) > E(S + 1), 
and the dependence on S is essentially linear with a small 
curvature. 

We conclude this section by discussing the different 
symmetry properties of even and odd chains. The mir- 
ror symmetry about a central site for odd chains can 
support a Neel-type antiferromagnetic configuration in 
the ground state [as indicated in Fig. ITUl' b). top picture]. 



For a number of bipartite molecular nanomagnets with 
antiferromagnetic Heisenberg interactions such as regular 
even wheels, modified even wheels, and grid molecules, 
the lowest lying spectrum was observed to consist of ro- 
tational bands or sets of states whose energies increase 
as E(S,q) = 7}AS(S + 1) + e(q), with appropriate gap 
A and "dispersion relation" e(q), where q is a suitable 
index or quantum numberi 28 ' 29 i 31 The lowest band (L 
band) and the set of higher-lying bands (E band) are 
distinguished according to a selection rule and the dif- 
ferent nature of the excitations associated to them. This 
structure of the excitations is reminiscent to that shown 
in Fig.CQfor the AFHCs, and the E(S) oc 5(5 + 1) energy 
dependence is indeed well obeyed for the odd chains, but 
as was demonstrated in the previous section pronounced 
deviations occur in the antiferromagnetic region for the 
even chains. In order to examine the even-odd effect the 
L&_E-band picture is therefore reanalyzed. 

In the Lh,E band picture the L band is described by 
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an effective Hamiltonian Hab of two collective spinsj^I 

Hab = AabSa • S B , (8) 

where Sa = J2ieA s * ano - &b = SieB s * are the total 
spin operators of the two sublattices A = {1,3,...} and 
B = {2,4,...}. The value of Aab can be determined 
using symmetry or mean-field arguments , 30 i 31 i 52 i 53 and 
is given for chains and rings in Eq. ([3]). The eigenvalues 
of Hab are easily determined to be 

E(S) = ^{S(S+1)-Sa(Sa + 1)-S b (Sb + 1)}, (9) 

where \Sa — Sb\ < S < (Sa + Sb)- For given values 
of Sa and Sb the energy spectrum consists of rotational 
bands. In the L band Sa and Sb assume their maximal 
values, and the slope is calculated to g(S) — 1. The 
next higher- lying states of Hab , which involve the states 
where either Sa or Sb is reduced by one, are related to 
the E band. However, they do not reproduce the energies 
of H because the spatial fluctuations, which give rise to 
dispersion of these excitations, are not accounted for in 
the Hab model [that is, e(q) is not obtained correctly] i^L 

Based on the observation that the energies follow an 
S(S + 1) dependence and that the calculated eigenval- 
ues and matrix elements become quantitatively exact for 
large s, the Hab model has been called classical or semi- 
classical. However, this notion is misleading as Hab re- 
tains a quantized spectrum. The basic underlying as- 
sumption of the model is in fact that the correlation 
length is infinite. The Hab model could then be re- 
garded as a symmetrized version of the bipartite Heisen- 
berg Hamiltonian, and in fact it is part of the greater 
class of the symmetrized effective Hamiltonians which we 
introduce next. 

With a unitary symmetry operator P, a translation 
operator for example, the symmetrized effective Hamil- 
tonians are constructed as follows: if a Hamiltonian H 
remains unchanged after j transformations with P, 

p- j HP j = H, (10) 

for example if P 3 — 1, then an effective Hamiltonian H c r 
is formed by adding up all j transformations of H, 

1 tl 

H eS =~Y / P- i HP i (11) 

This effective Hamiltonian commutes with P, since 
P~ 1 H c ffP = H c ff. A common eigenbasis can thus be 
found where the eigenvalues of the unitary P have mag- 
nitude 1, and one derives from Eq. (|11[) 

<*» ff | H cS | = (*™ ff | H | *Sr> , (12) 

since P % ^"ff) = il l^cff) f° r an y *• ^ then perturbation 
theory is applied to a non-degenerate eigenstate 
using H = H c ff + V , H c g is found to be exact in first 



order according to Eq. (|T2"j) . This establishes also an up- 
per bound for the ground-state energy of H, in any spin 
sector S if H is SU(2) invariant, as the higher orders in 
perturbation theory will result in a negative contribution 
by the variational principle. 

Hab is obtained if P is a sublattice permutation oper- 
ator, and it therefore follows that Hab is exact in first or- 
der for non-degenerate energy eigenstates or SU(2) mul- 
tiplets, respectively, and the L band in particular [see 
Figs[TJa), (b)]. For odd chains and even rings, it suffices 
to use a translation operator on one (the bigger) sublat- 
tice. j is then the size of the sublattice, and Aab = 
For even chains the two sublattices have to interchange. 

At this point it may be surprising that odd chains can 
be described by the Hab model, as the two collective 
spins Sa and Sb are necessarily of different size and 
therefore the problem becomes less symmetric. However, 
much more important is the fact that the collective spins 
obey the mirror symmetry for odd N, which is not the 
case for even N (where Sa and Sb are interchanged by 
this symmetry operation). For even chains on the other 
hand the Hab model would be expected to work well, 
in contrast to our findings. In a study on the N = 8, 
s = 3/2 chain it was argued that the L- and _B-band 
states mix producing deviations in the energies, which 
is forbidden in rings due to translational symmetry, at 
least in first order ^ The argument is obviously correct, 
but does not provide further insight, in particular on the 
even-odd effect. 

Except for energies with S < 2s in the even chains, 
Fig. [5] could suggest that the Hab model works other- 
wise well for chains. However, the slopes which approx- 
imate the L band in the ferromagnetic region in Fig. [2] 
are significantly larger than what is predicted by Eq. . 
This is not easily accounted for within the Hab model. 
The upper bound given by Hab for the energies in the 
L band is hence not very tight for both even and odd 
chains [see Figs.QJa) and (b)], giving room to convex de- 
viations from a S(S + 1) dependence, which in the case 
of the even chains indeed occur for S < 2s. 

It is also interesting to inspect the lowest energies in 
odd chains for S < S g and compare to those of Hab- 
Hab predicts a linear dependence on S, i.e., up to a 
constant E(S) = -±A AB (N- l)sS or S E (S) = S g -Sfor 
S < S g . A linear dependence is indeed observed in this 
spin sector as shown in Fig. |4l but the slopes are much 
lower than what is predicted by Hab ■ This is attributed 
to the neglected spatial fluctuations or dispersion of the 
excitations, which generally lowers the minimal energy. 

For even rings and other bipartite clusters it is well 
established that the Hab model becomes more valid the 
smaller N and/or the larger s is j 31 ' 54 Regarding chains, 
a similar tendency should be expected, which is not ob- 
served; e.g. the even-odd effect persists for large s as dis- 
cussed in Secini In contrast to rings the open boundaries 
in chains induce inhomogeneities, which are apparently 
not properly reflected in a model with infinite correlation 
length. 



7 



Taken together, the Hab model appears appealing at 
first sight based on its 5(5 + 1) energy dependencies, 
which are consistent with the energy spectrum in the 
AFHCs up to the even-odd effect displayed in Fig. [5J 
and because it is exact in first-order perturbation theory. 
However, the detailed analysis revealed its deficiencies, 
which appear to be connected to its neglect of spatial 
fluctuations or the assumed infinite correlation length. 



IV. PATH INTEGRAL REPRESENTATION: 
THE 0(3) NONLINEAR u-MODEL 

The 0(3) nonlinear er-model (NLSM) has the poten- 
tial to be suitable for short chains as it is based on the 
large spin, semiclassical limit, and can take inhomoge- 
neous states into account. In general, boundaries lead to 
residual topological contributions in the effective action 
for both half-integer and integer spin chains,— a concept 
that goes beyond the physics the Hab model is capable 
of describing. In the following we relate the concepts of 
the 0(3) NLSM with the numerical findings for the spec- 
tra of the antiferromagnetic Heisenberg model on short 
chains. In the limit of vanishing spatial fluctuations we 
indeed identify rotational bands in the spectra. For long 
chains the topological boundary terms have been inter- 
preted as remaining free spins at the edges of the chain 
(edge spins) . 43 i 44 We show that this interpretation is not 
supported by our numerical data and cannot explain the 
even-odd effect in the spectra of short AFHCs. On the 
contrary, we find that in short chains the alternating 
magnetization can no longer be separated from the uni- 
form magnetization, as it is assumed in many theories. 



A. Derivation of the 0(3) NLSM 

The 0(3) NLSM originates in an effective path integral 
representation of the AFHC in the low energy continuum 
limi t 12 i 13 (see also Appendix [C]) . The Euclidean action 
of a single isolated spin Sj with Hamiltonian H is best 
expressed in terms of spin coherent states reading 2 



that \£1{t + 8t)) — |fi(r)) is of order 0(5t) is used. How- 
ever, there is no physical reason why the difference of two 
paths at adjacent time steps should be a small quantity. 57 
Only for large s the overlap of two spin coherent states 
becomes infinitely small 2 , and therefore only in this limit 
the action in Eq. (IT31 becomes exact. 

For finite chains of length N the action A, Eq. dE 
generalizes to 



r P N 
A[0] = / d,T(n\H\fl) -isVw[fii] 



(15) 



where (Q\H\n) = Js 2 J2f =1 ^ ' ^i+i and = 
| ill, f22, . . . , Hiv) is a product of single spin states on 
each lattice site. In order to derive an effective field the- 
ory some assumptions must be made which are justified 
in the low energy and large s limit. Since it is expected 
that at least for short range order the chain exhibits siz- 
able antiferromagnetic correlations, the local spin field is 
decomposed into a Neel field n and a transverse canting 
field 1 



a(r) = (-ir +i n(x. i ,r)Wl 



+ -l(x„r),(16) 



where Xi = ia and a is the lattice constant. The fields 
n and 1 are assumed to be slowly varying and chosen to 
fulfill l(xj,r) • n(xj,r) = and |n(a^,-r)| = 1. Field 1 
contains the Fourier modes near and field n those near 
ko, the ordering wave vector. Thus 1 roughly represents 
the net magnetization, which is assumed to be small. 

Then Eq. ([TBJ is inserted in Eq. (JTSJ), L = Na is set 
and the continuum limit is taken. After some technical 
calculations (see Appendix |C|) the 0(3) NLSM is gener- 
ated: 



A 



L ,/3 

dx / dr 
I) Jo 



with v — 2Jsa, 7 = 2/s and 



.4 



top 



N 



Atop (17) 



(18) 



A\n] 



dT(n\H\ft) - isw[Sl], 



(13) 



where (3 = l/(ksT). States are labeled by a three- 
dimensional unit vector ft 2 = 1 defined by the eigenstate 
equation 17-Sj|f2) = s\fl). The kinetic energy of the spin 
enters in the last term in Eq. (|13[) . defined as 



isw[J7] = —is / drdu ft ■ (d T fl x d u fl), 
Jd 



(14) 



where D is the region bounded by the closed curve fi(r) 
on the unit sphere. The auxiliary variables u and r 
parametrize D. The action is unambiguously defined 
provided that 2s is an integer^ It should be pointed out 
that in deriving the action in Eq. (1131) the assumption 



In Eq. (fTT|) the field 1 has been integrated out giving 
1 = 4j(n x 9 T n). In order to understand the topological 
meaning of Atop, for even N Eq. (ITBl is rewritten as 



N/2 



A 



top 



is}^u[n(x2i-i)] -a;[n(a;ai)] (19) 



(=i 

«r f dx^- |n=n -d x ii. 



(20) 



After evaluating the integral in Eq. (|2"0"|) and taking into 
account an additional uncompensated Berry phase for N 
oddj££- one generally obtains 



A 



Lop 



i°- {4ttQ + w [n(L)] - (-l)^[n(0)]} , (21 
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where Q is integer valued and simply counts the winding 
number of the field n. For periodic boundary conditions 
and N even, Eq. (f2~Tj) reduces to A top — iOQ, the so- 
called (9-terixi ) 12 ' 13 where 9 = 2irs. A top depends only on 
the topology of the path. 

In general, with respect to the derivation of Eq. (IT71) for 
open boundaries, when performing the continuum limit 
for open chains additional boundary terms in the action 
are expected to be found. In some cases these terms can 
simply be included by introducing an effective length. 
For alternating sums though, the cases of summing over 
an odd or even number of terms need to be distinguished. 
Higher order terms in the bulk action are found which are 
assumed to be small [see Eq. (|C7I) in Appendix [C] . 

The topological term Atop = i8Q m Eq. (fTTj) is known 
to play a crucial role for the behavior in the thermody- 
namic limit. It enters the path integral with a factor 
e~ i27rs< 2, thus for integer spin it docs not affect the par- 
tition function of the chain, while for half-integer chains 
9 = 7r yields an alternating factor (— l)^. The renormal- 
ization group treatment of the action in Eq. (|17p with 
the constraint |n| = 1 leads to an increase of 7 which is 
at first independent of A top at large energy scales. With- 
out the topological term this leads to a strong coupling 
fixed point with a mass gap Ah ~ 0AJe~ ns for inte- 
ger spin chains, which has also been suggested from the 
calculation of the exact S matrix. 58 However, with the 
topological term for 9 = tt the renormalization behavior 
is altered and becomes difficult to analyze at small en- 
ergy scales, but it is known that half-integer spin chains 
are massless and belong to the universality class of the 
SU(2) Wess-Zumino-Witten model with topological cou- 
pling k = iSSriS in agreement with the Lieb-Schultz- 
Mattis theorem 6 ^ and the Bethe ansatz solution for the 
integrable spin-^ chain. This difference between integer 
and half-integer spin chains has become well-known in 
the literaturei 12 i 13 However, in the small chains that are 
considered in this work, the renormalization group flow 
is terminated by the finite length of the system at rather 
large energy scales, which exceed the mass gap, i.e., be- 
fore A t0 p makes a significant difference. Therefore, there 
is no need to distinguish between massive and massless 
theories in the following analysis. 



B. The Rotational Band in the 0(3) NLSM 

In order to connect with the Hab model, spatial fluc- 
tuations in the action are initially ignored. The problem 
then reduces to a (0 + I)-dimensional one and the action 
of Eq. (fTT)l taking into account Eq. (|2"Tj) is given by 

1 r 13 

A =o dT(d T n) 2 + ieaw[n], (22) 

2 Jo 

where / = N/AJ, e — for N even and e = I for N 
odd. For e = this yields the path integral of the 3- 



dimensional rigid rotator™, with energy levels 

E(S) = ^S(S + 1) (N even), (23) 

where S = 0, 1, . . . and each energy level is (2S + l)-fold 
degenerate. 

The case e = 1 can also be solved exactly, e.g., in the 
CP 1 representation by introducing an independent gauge 
field. 65 This auxiliary field can be set to zero, and the 
Lagrangian becomes the rigid rotator Lagrangian again. 
However, the variation of the gauge field generates a con- 
straint on the possible quantum numbers. One obtains 

E(S) = ±[S(S + l)-s 2 ] (iVodd) (24) 

with the allowed S now constrained to S = s, s + 1, . . . . 

For open boundaries there is an additional discrete lat- 
tice symmetry i — ¥ N+l—i, which is maintained by intro- 
ducing an effective length L = (N + l)a in the continuum 
model. The moment of inertia then reads / = (jV+l)/4J, 
and thus Eqs. (|2"31 and §Zty recover the slopes of the Hab 
model in Eq. ([3]) up to order 0(1/N 2 ). Taking into ac- 
count short range fluctuations in the Berry phase of a 
single spin^i or higher order operators in the effective 
action generally renormalizes the coupling constants. A 
similar mapping of magnetic molecular rings onto a rigid 
rotator model has been performed in Refs. (66ll67| . 

Neglecting spatial fluctuations in the 0(3) NLSM ap- 
parently renders its physics equivalent to the Hab model. 
However, if spatial fluctuations are included, the coupling 
to the rotational band has to be analyzed, which will al- 
ter Eqs. (|23|) and (|24|) . In particular, one would expect 
the spatial derivative term in Eq. (|17p to become more 
important once the total spin is excited, since the addi- 
tional net spin needs to be distributed along the chain. A 
discussion of this point with regards to the spin density 
will be given in Sec. IVIII However, as soon as the spa- 
tial derivative term is included in Eq. (|T7|) and the full 
Lagrangian is treated the spectra cannot be extracted as 
easily anymore, but it is possible to consider the simpli- 
fied case of edge excitations. 

C. Effective Edge-Spin Hamiltonian 

In Ref. [43| long s = 1 chains of even or odd length 
were modeled by the 0(3) NLSM coupled to s e d ge = \ 
edge spins. The constraint |n| 2 = 1 of the nonlinear a- 
field was relaxed by adding an artificial mass term and 
a repulsive An 4 interaction. On a mean field level the 
parameter A can be assumed to be small. Then the field 
n can be integrated out, and an effective Hamiltonian 
where only the edge spins couple to each other results^ 3 - 

Hedge = Jctt^'l ' S L J (25) 

where and s' L are spin operators representing the edge 
spins. Equation (|25[) is a valid approximation at en- 
ergies much smaller than the Haldane gap, for chains 
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much longer than the correlation length. The effective 
exchange interaction J c g ~ (— l) N e^ N ^ J between the 
edge spins is ferromagnetic for N odd and antiferromag- 
netic for N even, where £ is the spin-spin correlation 
length of the corresponding spin chain. The edge-spin 
picture hence gives rise to a pronounced even-odd dif- 
ference in the L-band energy spectrum at small values 
of S. However, this difference can only be derived for 
long chains, in contrast to the even-odd effect observed 
in Sec. [Hj which is a property of chains that are shorter 
than the correlation length. 

For long chains the existence of edge states is not re- 
stricted to integer spin. Based on Eq. (|2T1) and interpret- 
ing the residual Berry phase with free spins, edge spins 
of magnitude s e dge = § for s integer and s e d ge = o( s — f) 
for s half-integer have been proposed. 44 Half-integer spin 
chains have thereby been pictured as a continuum model 
of a spin- ^ chain coupled to two " impurity" spins of mag- 
nitude 3(5—5). Equation (|25p then predicts the following 
spectrum: for N even the ground state has S g = and 
the edge spins form a singlet. The lowest energy excita- 
tions can be constructed by exciting the two edge spins 
into an S = 1, 2, . . . , 2s e d ge state, and the excitation en- 
ergies are given by 



E(S) - E(0) = J eS S(S + 1) (N even). 



(26) 



For N odd, due to the ferromagnetic effective coupling 
Jeff < the lowest energy is obtained if the edge spins 
are coupled to their maximal value 2s e d ge (in case of half- 
integer spin the bulk spin-i chain contributes the addi- 
tional spin I to the total spin S g = s of the ground 
state). Excitations with S < S g are constructed by cou- 
pling the edge spins to lower spin values 2s e d ge — S' where 
S' = S g — S. The corresponding excitation energies are 

E(S g -S')-E(S g ) = \J cB \S , (S' + 1) (iVodd) (27) 

for S' = 1,2,..., 2s edge - For large system sizes the edge- 
spin picture has been numerically verified ! 15 ' 40 ' 43 ' 4 ^ 68 
However, for short system sizes it does not seem justi- 
fied to regard the edge states decoupled from the bulk 
states, as explained above. 

The coupling of spatial fluctuations to the bulk states 
for shorter system sizes may be described by a more gen- 
eral effective edge-spin Hamiltonian with the following 
ansatz H' — |j + Hi where Hi = A„l • (sj + s^) + A s n • 
[sj — (—l) N s' L ], which is appropriate if the coupling to 
the edge spins is weak. In case the two parameters A u 
and A s are both sufficiently small the field 1 can be inte- 
grated out, which recovers the edge-state Hamiltonian in 
Eq. (|25|) . However, it is so far unclear how to treat the 
effective model H' for the general case of short chains in 
order to extract the modified spectrum, which remains a 
task for future research. 

In Rcfs. [T5IJ45I ] it has been conjectured on the ba- 
sis of numerical density matrix renormalization group 
(DMRG) data that the edge-state picture stays robust 
up to very small N, even when TV becomes smaller than 
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FIG. 5: (Color online) Slopes G e d g e for odd (green squares) 
and even (blue circles) N. Solid symbols are for chains with 
s = 1, and open ones with dashed lines for chains with s = 2. 
The red lines represent the slopes G e d ge = J c s for large even 
N resulting from Eq. ^ with J cft = 0.35e _iv/6 J for s = 1 
and J eff = 0.05e _JV/33 J for s = 2 according to Ref. QJ]. The 
solid line is for s — 1, and the dashed one for s = 2. Lines are 
guides to the eye. 



the correlation length. In our data, however, we do not 
see signatures of edge states. First, for odd chains the 
lowest excitations E(S) for S < S g scale rather linearly 
with S in contrast to the quadratic behavior predicted in 
Eq. (f2"T)l (see Fig. 0] and the corresponding discussion in 
Sec. |H|. Interestingly, a linear dependence results from 
the Hab model. Second, for even chains the excitations 
at small S generally do not obey E(S) oc S(S + I) which 
contradicts the edge-spin prediction of Eq. (|26[) . In par- 
ticular, the deviation from quadratic behavior is observed 
up to ~2s, while edge states would correspond to the low- 
est s excitations [see Fig. Ufa)]. The edge states cannot 
be distinguished from the higher excitations up to ^2s. 
Finally, the scaling with chain length N of the first ex- 
cited state above the ground state does not follow the 
exponential decay of the edge-spin picture. To illustrate 
this let us define the corresponding (unnormalized) slopes 
in the E vs S(S + 1) diagram: 



A 



AB 



G 



edge 



. 2 



9(s) 



5(1) 



(N odd), 
(N even), 



(28) 



where g(S) is defined in Eq. ([2]). For N even, G e d ge 
is related to g AF by G edge = g AF A AB /2 (see also 
Fig. [3]), and for N odd, it is related to 5e by G e d ge — 
(N ~ l)A AB S E (s - l)/4 (see also Fig. g]). Fig. shows 
G e d g e for s — I and 2 for even and odd N. According 
to Eqs. (l2"rJl) and (|27l) . for long chains these slopes are 
given by G edge = Joff for N even and G edge = | J e ff|/s for 
N odd. The fit curves to the data for long even chains 
which thus have the form G e d ge — 0.35e _JV / 6 J for s = 1 
and G e d ge = 0.05e" Ar / 33 J for s = 2 where the fit pa- 
rameters are taken from Ref. [l5[ are also plotted. The 
deviations for smaller N are obvious. Furthermore, the 
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deviations appear to increase with increasing s, demon- 
strating that the edge-spin picture becomes less appro- 
priate with increasing s, in contrast to the observations 
for the even-odd effect of Sec. |TT1 which is present even 
for large s. 

In conclusion, the standard edge-spin picture cannot 
account for the even-odd effect. Including couplings of 
edge spins to both 1 and n in the NLSM would result in 
a more complete model that may remedy the situation. 
However, this would imply that uniform and alternat- 
ing magnetization become strongly coupled and can no 
longer be treated separately. 



V. COMPARISON TO THE SPIN 1/2 CHAIN 

As both the Hab model and the NLSM fail to account 
for the even-odd effect quantitatively, it may be instruc- 
tive to turn to the special case of s = 1/2 chains to 
examine if quantum effects play an important role. The 
spectrum of finite s = 1/2 chains is quantitatively very 
well understood, not only from the Bethe ansatz^ but 
also in terms of effective bosonic quantum numbers from 
bosonization , 51 ' 70 which establishes the s = 1/2 chain as 
an excellent reference. 

The description of the spectrum of the s = 1/2 chain 
in terms of bosonic quantum number o 51 ' 69 i 70 results in an 
almost equally spaced energy spectrum in the form of a 
conformal tower. There are corrections of order 1/N 2 and 
1 / (N In N) to the spectrum, but this effective description 
works well for N > 10. An L band can also be observed, 
except that in this case the lowest lying energy states of 
a given 5 are created by adding bosonic particles with 
zero momentum, and the number of bosonic particles is 
given by the S z quantum number, the projection of 5 
along the z axis. For s = 1/2 chains with N both even 
and odd the excitation energies in the L band are given 
by£i 



(a) 



E(S) 



N + l 



-S 2 



(29) 



up to higher order corrections in 1/N 2 and l/(N\n.N), 
where v — ttJ/2. Hence, there is no even-odd effect to 
lowest order in the excitation spectrum. There is a con- 
tribution to the ground state energy of order 1/N, which 
is positive for odd TV and negative for even JVjil but this 
is not related to the even-odd effect in Sec.|H] It is impor- 
tant to notice that in the s = 1/2 chain the E(S) depen- 
dence is predicted to be changed from the 5(5+1) behav- 
ior to a simple 5 2 behavior, analogously to the charging 
energy of a capacitor, which shows a quadratic energy 
dependence in the charge Q 2 . 

Interestingly, such an 5 2 behavior seems to agree bet- 
ter with the L-band energies for even chains with larger 
s, as shown in Fig. [6fa) for the case of N = 6 and s = 3/2 
corresponding to the Cr5 molecule. The S 2 behavior is 
consistent with the entire L band with -^qfj = 0.32 J or 
v ~ 0.71J according to the fit in the figure. In order to 
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FIG. 6: (Color online) (a) Energy spectrum of the N = 6, 
S = 3/2 chain as function of S(S + 1). The (red) line is the 
fit E(S) = 0.32S 2 . (b) Slopes g(S) x 25/(25-1) as functions 
of S for the N = 6 chain for different spins s (symbols). The 
dashed (red) line is the prediction of the Hab model. 



test this further, the slope of the L band in an energy 
vs. 5 2 diagram, which is given by g(S) x 25/(25 — 1), 
is plotted in Fig. [6jb) for N — 6 and different s. The 
5 2 dependence doesn't account for the even-odd effect, 
which in Fig. (Hb) is observed as the pronounced "dip" 
at small 5, but provides a significantly better average fit 
through the spectra compared to the 5(5 + 1) behavior 
in Fig. Ufa). However, this "success" of an 5 2 behav- 
ior does not explain why such an approach fails for odd 
chains, which obviously remain to be well described by 
the 5(5 + 1) behavior as shown in Fig. |2jb), nor does it 
give an independent estimate for j^rj- 

For the s = 1/2 chain it is also possible to pre- 
dict local expectation values, such as the alternating 
magnetizatio n 18 ' 71 i 72 and the dimerization along the 
chain. For example the alternating spin expectation val- 
ues in the z direction can be calculated for the highest 
weight states in the L bandit 



<*?) « (-1)" 



_ t siii(2Tr S z Xl /L) 
\J~L sm(irxi/ 1 L) 



(30) 



where L = (N + l)a is the effective length. Possible 
multiplicative corrections 7 -^ of order 1/ In L and higher 
order terms have been neglected here. The alternating 
order always decreases oc \fi near the edges^ The cal- 
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culation also implies an even-odd effect in the density: 
For the ground state of odd N chains with S z = 1/2 
there is a maximum oc \j\f~L in the middle of the chain, 
while the alternating order is zero for even N chains with 
S z = 0. The result shows explicitly that inhomogeneities 
are present and important over the entire chain, which 
were of course neglected in the Hab model. Spin densi- 
ties for higher s cases will be discussed in Sec. IVIII 

A similar calculation yields the alternating part of the 
nearest-neighbor correlation for states in the L band, 
which is dimerized 



> oc (-1) 



i cos(2n Sxi/L) 
L sm(jrXi/ L) 



(31) 



The dimerization becomes very strong and length inde- 
pendent at the edges (si • s 2 ) ~ —0.65 , but remarkably 
there is no pronounced even-odd difference, since the co- 
sine function near the boundary is independent of S being 
integer or half-integer. Corrections to the edge dimeriza- 
tion are small down to very short s = 1/2 chains and the 
correlation of the first two sites (si -S2) is much enhanced 
compared to the bulk value of —0.4431 both for even and 
odd N, despite the fact that only even chains could po- 
tentially lock into a dimerized ground state, while odd 
chains naively should not be able to support such a va- 
lence bond state. In fact, the difference of the expectation 
value (si • s 2 ) is only about 15% between chains of N = 4 
and N = 5. 

In conclusion, the analysis in this section shows that 
the even-odd effect described in Sec. QT] is not present 
in the quantum theory of the s = 1/2 chain. An even- 
odd difference in the spin density can be observed^ but 
this is related to symmetry properties, as also discussed 
later in Sec. IVIII However, strong inhomogeneities are 
observed in the chain and quantum effects cause the L 
band to be better described on average by a "charging 
energy" of the form E(S) oc S 2 . 



and nearest-neighbor correlation functions of the classi- 
cal AFHC, and compare them to the quantum results, as 
functions of the normalized squared total spin 



S ,, 



S(S + 1) 
5/(5/ + 1) 



(32) 



or S n = S/Sf in the classical case, where 5/ = Ns. 

The classical analog of Eq. ([T]) is constructed by intro- 
ducing unit vectors e.; = Sj/s, whose components com- 
mute in the limit s — > oo^^ The classical vectors can 
then be parameterized in spherical coordinates as = 
(cos 4>i sin#i, sinc/>i sin#i, cos 9i). Substitution in Eq. ([T]) 
minimizes the energy when 6i — O^+i-i and the nearest- 
neighbor relative azimuthal angles 4>i+\ — <f>i = 7T ) 77 i 78 
thus the spin configurations are planar as expected. The 
classical Hamiltonian then reads 



H = Js 



N-l 
2 \ ^ , 



+ 0i+i), 



(33) 



where 6*; <E [0,7r] for all i. Minimization of Eq. (|33|) 
gives the absolute ground state^Zr— Employing rota- 
tional symmetry, the lowest energy for arbitrary total 
magnetization S = s can be calculated by adding 

an external magnetic field term Hb = —BsY], cosdj in 
Eq. (f3"3"| , where B is directed along the z axis (the field is 
measured in units of 2/j.b in this section). The direction 
of the magnetization coincides then with the direction 
of the field and S — sJ2i C0S ^i- By tuning B the zero- 
field energies E(S) can be calculated for all values of S 
by subtracting the magnetic energy at the end. For odd 
chains, configurations with magnetization less than the 
value of the absolute ground state S < S g = s are not 
accessible this way, as E(S) decreases as function of 5 in 
this regime. The calculation of these states is performed 
by adding a term Hk = KS 2 with K > 0, which favors 
states with minimal S. 



VI. CLASSICAL ANTIFERROMAGNETIC 
HEISENBERG CHAIN 

The classical AFHC, where the spins in Eq. (JlJ arc 
treated as classical objects (vectors), is a good approx- 
imation of the quantum model for large s. It ignores 
quantum fluctuations but fully retains spatial fluctua- 
tions, and thus allows to study the importance of the 
latter. This model has already been studied within the 
context of artificial nanostructureS ) 73 ' 74 and it has been 
found that even chains always have a coplanar and non- 
collinear ground state. For odd chains the situation is the 
same, except that for magnetic fields below a critical field 
the lowest energy configuration is ferrimagnetic, and an 
analytical expression for the critical field was founds 
This difference reflects the different total ground-state 
spins S g = and S g = s for even and odd chains. In 
this section we present the lowest energy, spin density 



A. Slopes and Energies 

The slope in the classical case is determined as g(S) = 
-^dE/d{S 2 ) = ^sdE/dS, and is related to the in- 
verse magnetic susceptibility (S/B)" 1 : According to the 
Legendre transformation E(S) = Eb + BS, the magne- 
tization S is given as S(B) — —dEs/dB (at T = 0) and 
the field as B(S) = dE(S)/dS. One thus finds 

B(S) = A AB Sg(S). (34) 

The inverse of this equation gives the magnetization as 
a function of field, S(B), which implies that AABg(S) is 
the reciprocal susceptibility (S/B)^ 1 (see also Sec. HI)) . 

The L-band energies of chains with lengths N = 6 
to 1 1 are displayed in Fig. [TJa) as functions of S 2 , and 
the corresponding slopes are presented in Fig. EJb). For 
even chains, the classical slope g is small at small S and 
increases rapidly with increasing S (which comes about 
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FIG. 7: (Color online) (a) L-band energies E(S) of the clas- 
sical AFHC scaled with J(N — l)s 2 for chains with lengths 
TV = 6 to 11 for smaller values of S„. For the odd chains 
the energies with S < s do not belong to the L band but to 
configurations with magnetization less than the one of the ab- 
solute ground state, and are included here for completeness. 
The arrow points towards the "kinks" in the energies at the 
fields where S = s. The inset shows the same figure for the 
whole S„ range, (b) The slopes g of the L band of the classi- 
cal AFHC as function of S„ for chain lengths TV ranging from 
6 to 11. 



effectively like an increase of the external magnetic field 
B and can be thought of in these terms). In contrast, 
the slope (for S > S g ) for the odd chains starts off at 
a higher value compared to the even chains, and shows 
a significantly weaker dependence on the total magne- 
tization. The slopes for both the even and odd chains 
become comparable and weakly varying for S% « 4/TV 2 
or S Ri 2s. In Fig. [5] the energy spectra E(S) of the 
quantum AFHC, scaled with the energy of the ferromag- 
netic state, J(TV — l)s 2 , are shown for the TV = 6 and 7 
chains for s ranging from 1/2 to 7/2, and the classical re- 
sults are also shown. The corresponding slopes g(S) are 
presented in Figs. [9ja), (b). The quantum energies and 
corresponding slopes approach the classical ones with in- 
creasing s. Although convergence is relatively slow, in 
both the even and odd chains the classical and quantum 
slopes exhibit very similar features. Most importantly, 
the strong down-bending in the slope for the even chain 
at small values of S 12 (or S), which is the hallmark of 
the even-odd effect of Sec. UH is also present in the clas- 



FIG. 8: (Color online) L-band energies E(S) scaled with 
J(TV - l)s 2 for the (a) N = 6 and (b) N = 7 chain, for s 
ranging from 1/2 to oo (symbols). Results are shown as func- 
tion of S 2 . For TV = 7 (odd chains) the energies with S < s do 
not belong to the L band but to configurations with magne- 
tization less than the absolute ground state, and are included 
here for completeness. 



sical system. This implies that spatial inhomogeneities 
must be the leading mechanism of the even-odd effect, 
while quantum fluctuations give quantitative corrections 
to this phenomenon. 



B. Local Magnetization and Correlations 

To study the spatial inhomogeneities further local 
quantities are considered. The local magnetizations of 
the spins along the direction of S (with S > 0) or spin 
densities e\ = cos 9i and the nearest-neighbor correlation 
functions ei-e i+ i are shown in Figs. EJc), (e) and^d), (f) 
for TV = 6 and 7 respectively. The related spin configu- 
rations are presented in Fig. [TU] for different values of S 2 . 
e\ and e,--ej-|_i obey the mirror symmetry of the chain. In 
even chains the local magnetization ef is zero everywhere 
in the ground state as the spins align perpendicular to S. 
For very small S (or magnetic field B) the ef are almost 
perpendicular to S, to optimally preserve their exchange 
energy. With increasing S the outer spins ef and e z N have 
the largest projection on S among all spins and gain the 
most magnetic energy, while their nearest neighbors e| 
and e z N _ 1 turn against the field. This configuration al- 
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FIG. 9: (Color online) Behavior of the classical AFHC for N = 6 (top) and N = 7 (bottom) as function of (a,b) The 
classical slopes where s — > oo (red lines), along with the slopes of the quantum AFHC (symbols) of Fig. [2] are shown. (c,d) 
Spin density ef = cos 8i along the direction of total spin S. The numbering of the spins start at the edge of the chain. The 
spin density is mirror-symmetric with respect to the center of the chain. (e,f) Nearest-neighbor correlation functions e» ■ ej+i. 
They are mirror symmetric with respect to the center of the chain. 



lows a net magnetization at low exchange energy cost, 
and the edge spins are in fact very quickly magnetized, 
which implies a large magnetic susceptibility or a small 
slope g. The situation is quite different for odd chains 
since the ground state is in a ferrimagnetic configuration 
and the outer spins are already fully aligned with the to- 
tal magnetization S = s. In order to magnetize the chain 
further the spins on the odd sites decrease their local 
magnetization ef , and this allows the spins on the even 
sites to increase their magnetization. This magnetization 
process costs more energy and is less efficient than in the 
even chain. Hence the susceptibility is smaller and the 
slope g(S) is larger for odd chains. This is also reflected 
in the local correlations in Figs.[^e), (f) where the corre- 
lation ei ■ e2 of the first bond increases more strongly for 
odd chains than for even chains, which in turn requires 
more energy. The markedly different slopes at small S 
in the even and odd chains are hence related to the high 
susceptibility of the outer spins in the even chains to- 
wards magnetic fields, providing an intuitive picture of 
the even-odd effect. 

Remarkably, at S„ ~ A/N 2 or S w 2s the differences in 
correlations and local magnetizations between even and 
odd chains start to disappear, and the slopes g(S) in 
Fig-IHa) and (b) become comparable [see also Fig.rTJb)]. 
This crossover region S ~ 2s is marked by a vertical 
dotted line in Figs. [HIc),(d) and the corresponding spin 
configurations are depicted in Fig. I10I For S larger than 
2s the interior spins exhibit nearly identical ef for both 



even and odd chains, while the outer spins have larger 
local magnetization. 

C. Analytical Results 

Further insight is provided by the analytic calculation 
of the local magnetization along the chain in the limits 
NB — > and NB — > oo, respectively. Using a small 
angle expansion of Eq. p3[) for infinitesimal fields B we 
arrive at a set of coupled equations, which can be solved 
analytically for the local magnetization of classical chains 
with arbitrary even N 

ef = 47; [l + ( " 1)4(2l ~ 1 ~ Ar)] (35) 

for i < N/2 (e z N+1 _ i = ef due to mirror symmetry). 
The local magnetization hence decays linearly from the 
edges into the chain. The total magnetization is obtained 
as S — NB/(2J), where it should be noted that both 
the uniform and alternating parts in Eq. Q35[) contribute 
equally to it. This is remarkable, since it implies that 
the alternating part due to the boundary condition is 
affecting a thermodynamic quantity, i.e., the edge effect is 
of order N and cannot be extracted from standard finite- 
size scaling. The energy is E(S) = -(JV-1) J+(J/N)S 2 . 
For the slope thus holds 
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This is almost a factor 2 smaller than the prediction 
g = 1 of the Hab model or the slope gp in the ferro- 
magnetic region S n — > 1 [Eq. ([SJ], and thus explains the 
strong reduction of g^S 1 ) for even chains analytically (but 
does not explain a crossover at S 2s). For odd chains 
we could not find a closed analytical solution of the lin- 
earized equations. 

The local magnetization can also be calculated approx- 
imately from an effective hydrodynamic theory in a semi- 
infinite chain in finite fields (i.e. NB ^> Js)^ - which can 
also be derived from the classical version of the NLSM in 
Sec. HVl This results in the following expression for the 
local magnetization 80 



^-(-lyexpH/^) 



dq cos(jg) 
^ 2^ U 2 s 2 q 2 + B 2 



(37) 



without any adjustable parameters (i < N/2). Here, we 
introduced the quantity 



2Js 
~B~' 



(38) 



which defines a characteristic length in units of the lat- 
tice spacing. The local magnetization decays exponen- 
tially into the chain with a length scale £b, which de- 
pends on the field B. Interestingly, the prefactor of the 
alternating part is unity and independent of J and B. 
The field therefore does not determine the strength of 
the alternating response but only its range £b- In this 
case, the edge effect is also large but not of order TV, 
and the thermodynamic contribution can be extracted 
using finite-size scaling, in contrast to Eq. (|35|). The lo- 
cal magnetization of a N = 100 chain at intermediate 
field is shown in Fig. 111! and good quantitative agree- 
ment is found with the exact classical result [for very 
short chains Eq. (f3"T|) holds only qualitatively, see below] . 
In first order in (NB) -1 the total magnetization is ob- 
tained as S — NB/(AJ), and £b thus varies with S as 
£b oc S^ 1 . Since S < Sf (magnetic field < saturation 
field) the limit £b > 1/2 is implied. For relatively large 
S, when £b <C N, the alternating part becomes located 
near the edges, and the local magnetization becomes es- 
sentially homogeneous in the interior of the chain, as it 
is also observed qualitatively in the spin configurations 
shown in Fig. [TU] for S 2 — 0.3. The slope is determined 
as g(S) — 1, as in the Hab model, and consistent with 
the exact ferromagnetic slope gF which is approached for 
TV — > oo . This finding may serve also as a measure of the 
quantitative accuracy of Eq. (|37l) . 

Eqs. (|35|) and (|37|) describe two completely different 
physical regimes, and combined they provide a qual- 
itative description of the crossover. At small fields 
B <C Js/N the characteristic length £b exceeds the chain 
length and there is full interference between the edges, 
while at large fields B ^> Js/N the length £g is much 
shorter than N and the two edges act independently. Ac- 
cording to Fig.|H]the crossover occurs when S ~ 2s, which 
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different values of extracted from Fig. [9] (c) and (d). 
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FIG. 11: (Color online) The local magnetization ef along an 
even classical chain of length N = 100 for a field of B = 
0.2J (solid symbols) compared to the theoretical prediction 
in Eq. ()37p (open symbols). 



translates into N m 4£g [which also implies N 3> 2 for 
Eq. (|3"T|) to describe quantitatively the large-field regime, 
which is not always fulfilled in very short chains]. 

In the thermodynamic limit N — > oo at very low fields 
B — » both Eqs. (|35p and (|3"7} appear to be valid, but 
give contradictory results. This discrepancy is resolved 
if one takes care of the order of limits. If NB — >• 
only Eq. (|3"5j) is applicable while Eq. (|3"?| only holds if 
NB — > oo. Therefore, the thermodynamic limit N — > oo 
and the zero-field limit B — > do not commute in the 
classical model with edges, which has also been observed 
for impurity effects in higher dimensions.— _ — 

In this section it has been shown that the difference 
in the i-band behavior between even and odd chains 
is captured by the classical AFHC model. The ana- 
lytic calculations for the spin densities naturally suggest 
a crossover in even chains at the onset of interference 
between the edges, which the numerical results show to 
occur at S ps 2s or N s» 4£g. Below this magnetization 
the edge spins of even chains can be magnetized with a 
low energy cost, leading to a reduction of g(S) predicted 
in Eq. (|36|) . In odd chains on the other hand the ferri- 
magnetic configuration at small S prevents an easy mag- 
netization of the edge spins, leading to a larger slope g(S) 
in the numerical results. It should be noted that the edge 
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FIG. 12: (Color online) Spin density (sf) and nearest- 
neighbor correlations (s; ■ s;+i) of the s = 1 AFHC (black 
solid circles), (a) Spin density in the lowest S = 1 multi- 
plet of the N = 16 even chain and (b) for the S = 1 ground 
state of the TV = 15 odd chain, (c) Correlations in the S = 
ground state for TV = 16 and (d) in the S = 1 ground state 
for N = 15. For comparison in panels (a) and (b) the spin 
density of the corresponding Hab model is shown (blue open 
circles). The nearest-neighbor correlations in the Hab model 
are equal to —1.125 (not shown). Lines are guides to the eye. 



spins here are classical and are hence not related to the 
quantum edge spins of the NLSM (Sec. IIV C[) . At mag- 
netization above 2s the alternating spin density localizes 
at the edges, and the distinction between even and odd 
chains largely disappears. 



VII. SPIN DENSITY AND CORRELATION 
FUNCTIONS 

After having shown in the previous section that the 
even-odd effect can be rationalized with the help of the 
classical spin densities and correlation functions, these 
quantities will be briefly examined for the quantum 
AFHC in the antiferromagnetic region. The spin den- 
sities and correlation functions have to be symmetric un- 
der the parity operation with respect to the center of the 
chain. This leads to an obvious difference in the spin den- 
sity or wavefunctions between even and odd chains. For 
even chains the parity operator interchanges the two sub- 
lattices, flipping the spins of each sublattice. The symme- 
try competes with the antiferromagnetic order and leads 
to having the same nearest neighbor (sf ) around the cen- 
ter, which is very small in magnitude. For odd N there 
is no such restriction. 

A specific example is shown in Fig. 1121 The spin den- 



sity of the s = 1 chains with N = 15 and 16 is plotted for 
the S = 1 lowest state (which is the ground and the first 
excited state respectively), showcasing the differences for 
even and odd N (note that for the ground state of the 
N = 16 chain (sf) = for all spins). The spin den- 
sity is weaker at the center, while it increases approxi- 
mately linearly going towards the edges, where spins are 
less bound. The predictions of the Hab model are also 
plotted in the two figures, and they miss the main fea- 
tures, even though they exhibit a difference between even 
and odd chains. For the odd chain the Hab prediction 
clearly shows the antiferromagnetic order, while for the 
even chain it is uniform (and hence small in magnitude) . 

For comparison, wc shortly comment on the situation 
in long AFHCs. In Ref. [43| the spin density in the lowest 
5 = 1 state was numerically calculated with the DMRG 
method for s = 1, N = 100. A very good agreement 
has been found with an exponential decay away from 
the chain ends, resulting from s e d ge = 1/2 edge spins. 
Additionally, data for higher spin states S > 1 confirmed 
the analytic picture of edge states and dilute boson-like 
bulk magnons in long chains^ However, for the s = 1 
chain the correlation length is about 6 sites, and the end 
spin wavefunctions protrude accordingly from each end 
of the chain. The data shown here are for N — 15 and 16 
which is about two times the correlation length, and the 
effective description by independent quantum edge spins 
is apparently not appropriate. 

Looking at the nearest-neighbor correlation functions 
for the ground states of the N—15 and 16 chains 
[Fig. I12f c).(d)]. only weak differences between even and 
odd chains are observed. Correlations are maximal at 
the edges, where the relatively loosely bound spins have 
more freedom to minimize correlations in their vicin- 
ity. The strength of the nearest-neighbor correlations 
decreases towards the center. The difference between 
the even and the odd chain is seen in the central re- 
gion, where the correlation oscillates in strength with 
position for A = 16, similarly to what happens for the 
spins further out [Fig. fT2Td)]. In contrast, for N = 15 
the strength of the central correlations does not oscil- 
late much with position [Fig. H2T c)]. The Hab model 
predicts for both cases uniform nearest-neighbor correla- 
tions equal to —1.125 and completely misses the central 
features. 



VIII. 



VALENCE BOND STATES 



The Hab model can describe the physics of the AFHC 
when Neel-type correlations prevail. Its main deficiency 
is that it doesn't account for spatial fluctuations, lead- 
ing to an infinite correlation length in this model and 
small nearest neighbor quantum entanglement. The 
valence bond solid (VBS) is a complementary descrip- 
tion, where strong (singlet) entanglement between near- 
est neighbors is built in, and correlations are exponen- 
tially decaying j 49 ' 50 i 84 In contrast to the Hab model the 
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FIG. 13: Sketches of valence bond states for N = 6 and 7 
chains. Each column represents a spin s, each circle a spin 
1/2, and each line a singlet bond, (a) Even and (b) odd 
completely dimerized s = 1/2 chain, (c) Even and (d) odd 
VBS state in an s = 1 chain, (e) Even and (f) odd partially 
dimerized s = 3/2 chain, (g) Even and (h) odd VBS state in 
an s — 2 chain. 



VBS states also explicitly contain spin degrees of free- 
dom near the edges, and might hence better approximate 
the spatial fluctuations relevant for the even-odd effect. 
We therefore compare and combine quantum VBS states 
with Neel-type Hab states in order to understand better 
which effect plays a more dominant role. 



A. Construction of Valence Bond States 

VBS states were originally introduced as translation- 
ally invariant ground states of exactly solvable integer 
spin models with an excitation gap , 49 ' 50 rigorously ex- 
emplifying the Haldane phas o 12 i 13 for the first time. In 
general, valence bond states are formed by replacing the 
spin s operators by symmetrized 2s spin-1/2 objects on 
each site, and then coupling pairs of spin-1/2 objects on 
different sites to form singlets^ Many different valence 
bond states can be constructed for a particular system in 
this way depending on the coupling scheme of the spins, 
forming an overcomplete basis of the Hilbert space in the 
singlet sector^ As a simplification, we here restrict the 
valence bonds to connect nearest-neighbors evenly to the 
right and the left as shown in Fig. [13] Because the spin- 
1/2 objects are symmetrized at each site, the resulting 
correlations remain extended over a correlation length of 
several spins. For integer spin chains this construction 
leads to a unique, and in case of periodic boundary con- 
ditions, translationally invariant state called a VBSj 49 i 50 
In case of half-integer s the number of singlets between 
nearest neighbor sites of the valence bond wavefunction is 
different for two successive pairs, see Fig. [TBTe.D. there- 
fore translational invariance is lost. The depicted state 
for each half-integer s case is complemented by a state 
where all bonds are shifted one lattice spacing to the 



right. More generally, a half-integer spin VBS can be 
regarded as an integer spin VBS with an additional spin- 
1/2 chain. However, according to the magnetization pro- 
file in Eq. (|3"U|) the residual free spin for odd N spin-1/2 
chains is located mostly in the center of the chain i.e. not 
as depicted in Fig. fToTb.f) near the edge. The approxi- 
mate ground state for half-integer s is therefore formed 
by all states where the residual free spin is delocalized. 
However, using a suitable parent Hamiltonian based on 
projection operators a unique trial VBS state can be de- 
fined as will be shown below. 

Using the original idea from Affleck, Kennedy, Lieb 
and Tasaki, parent Hamiltonians with nearest neighbor 
VBS wavefunction as the ground state can be constructed 
for integer s by projecting out all parts where neighbor- 
ing spins couple to a total spin less than s. In case of 
open boundaries, this results exactly in the ground states 
shown in Fig. fTBT c.d.g.h) with s unpaired spin 1/2 objects 
at each end. These edge spins can form total spin multi- 
plets ranging from to s, thus the ground state is (s + 1) 2 
times degenerate (including the degeneracy with respect 
to S z ). Parent Hamiltonians can also be constructed for 
half-integer spin chains^, as will be shown below for the 
case s = 3/2. 

In order to define parent Hamiltonians with exact VBS 
ground states it is useful to define (non-normalized) pro- 
jection operators acting on sites i and which project 
out all states with total spin (sj 



s m ) 2 < F: 
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The constant Kp > is conveniently fixed such that the 
prefactor of Sj • Sj+i equals 1. The bond projection oper- 
ators can be easily expressed as polynomials of Sj • Sj+i^ 
The parent Hamiltonians for s = 1 and 2 are ther* ^ 49 ' 50 
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For half-integer spin the situation is more complicated, 
because the total spin of two neighboring spins is obvi- 
ously alternating in the VBS states in Fig.fToTa.b.e.f) and 
cannot be fixed to a constant. One solution is to use al- 
ternating parent Hamiltonians. In particular for s = 3/2 
one can choose 



rrl2 

U s=3/2 



H 



s=3/2 



= E( p 



2i,2i+l 



P. 



2i+l,2i+2 



), 



(42) 



17 




2 4 6 8 10 12 14 

N 



FIG. 14: (Color online) Variational energies Evbs and Eab 
of the AFHC Hamiltonian for the VBS states sketched in 
Fig- EH (c)-(h) (green solid symbols) and the ground states of 
the Hab model (blue open symbols). The energies are shown 
with respect to the exact ground state energy E g of the AFHC 
for s = 1 (squares), s = 3/2 (circles), and s = 2 (stars), as a 
function of chain length N. Lines are guides to the eye. 



with two different ground state wavefunctions: | "£'12) 
with one singlet bond between the first two sites and 
|4 , 2i) with two singlet bonds between the first two sites. 
For even N only |^2i) is a reasonable VBS trial state, 
while for odd N both states are equivalent, so that a par- 
ity symmetric combination of the two must be formed. 

In the following we analyze the overlaps and expec- 
tation values of the corresponding VBS states \^vbs} 
depicted in Fig.[T3tc)-(h). The states \^vbs) can be nu- 
merically calculated as ground states of the parent Hamil- 
tonians using the iterative power method and a projec- 
tion onto the S subspace of interest or alternatively using 
an iterative method described in Appendix [B] 



B. Comparison of the Hab Model with the VBS 
Model 

In Fig. Q3] the ratio of the variational energies Eab = 
(y AB \H\y AB ) and E VBS = (^vbs\H\^vbs) over the 
exact ground state energy of the AFHC of Eq. ([1]) are 
plotted for s = 1, 3/2, and 2. I^ab) is the ground state 
of the Hab Hamiltonian of Eq. @, and \^>vbs) is the 
ground state of the corresponding VBS parent Hamilto- 
nians of Eqs. ([iO]). (|41~1) . and (|l2"|). The accuracy of the 
variational energy of the ground state of the Hab model 
drops off quickly with TV, however for very small N the 
Hab variational energy is better than the VBS varia- 
tional energy. Increasing s also improves the quality of 
the Hab variational energies, which agrees with the ex- 
pectation that the Hab model is best suited for small 
N and large s. It should be mentioned here that in the 
variational energy the even-odd effect appears to be re- 
versed: ground state energies for even N are on the aver- 
age slightly better approximated than for odd N chains 
by Hab, which is opposite to what would be expected 



from the behavior of the excited states as described in 
Sec. HU For the VBS energy ratios the variational energy 
for s = 1 starts out relatively poor, but then improves 
with increasing N. In contrast to the Hab model vari- 
ational energies, the energy is generally estimated well 
also for large N by \^vbs)- 

In Fig.[15]the overlaps of \^ab) and \^vbs) with the 
ground state of the AFHC of Eq. (JXJ) are shown. These 
overlaps largely confirm the picture discussed in the pre- 
vious paragraph. For small N and especially larger s 
the Hab model has a slight advantage over the VBS 
ground state, but then its overlap drops off quickly with 
N. Again the Hab overlaps are slightly better for even 
N than for odd N. The overlap of \^vbs) on the other 
hand is much less dependent on s and also drops off slower 
with N. This shows that local quantum entanglement is 
important for any s and N, while the Neel-type order of 
Hab is relevant for large s and small N. The overlap of 
the two ground states I^ab) and \^vbs) is also plotted 
in Fig. 1151 Interestingly, both models give very similar 
wavefunctions up to JV = 4, as can be concluded from 
the large overlap values. 

To improve on the quality of the variational approx- 
imation both the Hab and VBS models were simul- 
taneously taken into account, by forming a trial state 
as a linear combination of the two corresponding wave- 
functions, namely the optimal wavefunction |$ pt) — 
&\^ab) + b\^VBs) (the notation always implies normal- 
ization) . Its overlap with the AFHC ground state is also 
plotted in Fig. [15] for the optimal combination of the 
variational parameters a and b. The overlap improves 
in comparison with the two individual wavefunctions but 
still decreases with N, with a weak dependence on s. In 
the insets of Fig. [I5]the optimal ratio a/(a + b) is plotted. 
The overlap decrease with N shows the importance of the 
VBS state for longer chains. It generally increases with 
s, and the Hab model gains more weight as the increase 
of s makes the AFHC ground state more Neel ordered 
and less entangled. 

Comparing the approximation of the ground state en- 
ergy (Fig. rH| with the overlap of the ground state wave- 
function (Fig. [T5]) for the VBS model, the former hardly 
worsens with N, while the overlap of the wavefunctions 
decreases. This is due to the fact, that the overlap of the 
VBS model wavefunction with the AFHC wavefunctions 
of other low lying energy levels is still significant. Hence 
the VBS wavefunction mostly mixes with the low lying 
AFHC energy levels. 



IX. SUMMARY AND CONCLUSIONS 

In this work the structure of the lowest S excitations in 
AFHCs of relatively short length N but relatively large 
spin magnitude s has extensively been studied by con- 
trasting the results of a broad array of theoretical tools 
and approaches. The results of this paper are of relevance 
from at least three perspectives. 
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(a) s=1 — <*J*> 2 (b) s=3/2 (c) 




N N N 

FIG. 15: (Color online) Squared overlap of the ground state wavefunction of the AFHC with the VBS state ({9vbs\^} j 
green squares), with the Hab ground state ({^ ab\^} 2 , blue circles), and with an optimized linear combination (|$ pt) = 
ab) +b\ i &vBs}, black stars) for the spin values (a) s=l, (b) s=3/2, and (c) s=2 as function of chain length N. The relative 
overlap (^vbs\^ab) 2 is shown for comparison (red diamonds). In the insets a/(a + b) is plotted (black squares). The red line 
corresponds to a value of 0.5. Lines are guides to the eye. 



Quantum vs spatial fluctuations. First of all, the 
findings further our understanding of the physics in the 
AFHC model. It has been demonstrated that there is a 
distinctive even-odd effect in the dependence of the low- 
est energies E(S) in each total-spin sector on S or the L 
band in short chains. The effect is markedly different to 
the established even-odd effect in long chains, which is 
well understood in terms of the quantum edge-spin pic- 
ture; the arguments were given in Sec. II V Cl The different 
physics found in these two regimes justifies a distinction 
into short and long chains, which represents a major find- 
ing of this work. The described even-odd effect manifests 
itself in the antiferromagnetic region of the L-band spec- 
trum (low S), but not in the ferromagnetic part (high S). 
In the antiferromagnetic region even-odd effects can also 
be noticed e.g. in the ground-state spin, the spin density, 
and the nearest-neighbor correlation functions. These 
are however straightforwardly explained by the different 
symmetry properties of even and odd chains. The even- 
odd effect focused on in this work in contrast is not as 
trivially traced back to the different symmetry properties 
of even and odd chains. 

To elucidate the physics giving rise to this effect, differ- 
ent models were investigated, and the AFHC model was 
firstly compared to the Hab model. Phenomenologically, 
the Hab model appeared as a promising candidate since 
it naturally produces an E(S) oc S(S + 1) energy de- 
pendence and an L&i?-band structure, as approximately 
observed in short AFHCs (Fig. [1} . While for odd chains 
the Hab model describes the L band surprisingly well in 
the full range of S values, with a slight renormalization of 
the slope g(S) or the effective gap Aab, it fails to do so 
for even chains. The deviation is most pronounced in the 
antiferromagnetic region for S < 2s, which is the even- 
odd effect, but also in the ferromagnetic region the de- 
viation is significant. For even-membered antiferromag- 
netic Heisenberg rings the energies predicted by the Hab 
model were previously shown to become more accurate 
the larger s, and the Hab model was hence considered 



(semi-) classical in nature^ Our results on the AFHC 
correct this view and point to the fact that the main 
characteristics of the Hab model is the neglect of spatial 
fluctuations or implicit assumption of an infinite corre- 
lation length, which consistently explains our findings. 
For instance, the spin densities are more homogeneous 
in odd than even chains suggesting a better accuracy of 
the Hab predictions in the odd chains. Also, that the 
Hab model reproduces energies and transition matrix el- 
ements extremely well for rings is now expected from the 
fact that in rings the L-band states exhibit homogeneous 
spin densities by symmetry. These trends for rings, odd 
and even chains also manifest themselves in the ferro- 
magnetic region, as characterized by the slope gp. For 
rings one finds gp = 1, which coincides with the predic- 
tion of the Hab model (g = 1), while for chains (with 
N > 3) gp > 1 holds, with the larger discrepancy for the 
even chains. 

The predictions of the 0(3) NLSM for short AFHCs 
were also analyzed. Interestingly, with neglected spatial 
fluctuations the Hab model is reproduced, which under- 
pins the role of spatial fluctuations and establishes the 
theoretical basis of the Hab model. More importantly, 
the analysis demonstrated that the even-odd effect is not 
easily reconciled within the NLSM. It in fact showed that 
the usual assumption of describing Neel order and uni- 
form canting as separate, weakly coupled degrees of free- 
dom, which is exploited in or is even at the heart of many 
theories such as bosonization, hydrodynamic theories, or 
those based on the NLSM, fails in short chains. In par- 
ticular, it demonstrated that the even-odd effect is dis- 
tinct from the even-odd effects due to the quantum edge- 
spin model established for long chains. The latter was 
furthered by an analysis of the VBS wavefunctions for 
AFHC systems. As e.g. demonstrated by the analysis of 
the spin densities and correlations, differences between 
half-integer and integer spin chains are small for short 
chains. The VBS results thus suggest that in short chains 
(with s > 1/2) the integer-spin VBS part in the total 
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VBS wavefunction is more relevant to the physics than 
the additional half-integer spin VBS part present in half- 
integer spin chains. Somewhat surprisingly it was found 
that the VBS and Hab wavefunctions approximate the 
exact wavefunctions nearly equally well (or equally poor) 
for small chains with relatively large spin magnitudes s. 
The Hab and the quantum VBS models capture differ- 
ent aspects of the wavefunctions; each model has hence 
its strengths and weaknesses with no clear advantage for 
one over the other. 

Finally, the AFHC was also analyzed in the classical 
limit. The classical model describes the general trends 
and in particular the even-odd effect very nicely, as shown 
by the numerical results, leaving little doubt that it cap- 
tures the essential physics. As main result it demon- 
strates the importance of spatial fluctuations in the even- 
odd effect. Qualitatively, the even-odd effect can be re- 
lated to the spatial inhomogeneities introduced by the 
spins at the edges and their larger response to weak ap- 
plied magnetic fields in the case of even chains [i.e. a 
smaller slope g(S)j. At a quantitative level significant de- 
viations remain unexplained, which reflects the fact that 
for the considered spin magnitudes the classical limit is 
not yet reached and quantum fluctuations still play a 
significant role. As a striking, yet so far unexplained 
consequence, the L band of odd quantum chains is well 
described by the Hab model, which predicts the slope 
better than the classical model in Fig. HJ^b). 

Physical regimes in the AFHC model. Having es- 
tablished fundamentally different behavior for short and 
long chains, the question arises where the crossover be- 
tween these regimes is located. It was first argued by 
Haldane that in the thermodynamic limit s determines 
the physical behavior, leading to a gap Ajj ~ 0.4Je~ 7rs 
in the excitation spectrum for integer s, while half- 
integer spin chains have a linearly dispersing excitation 
spectrum^ - — In the framework of the renormalization 
group treatment of the NLSM in Eq. (TT71) both integer 
and half-integer spin chains in fact show the same in- 
crease of the dimensionless coupling constant 7 in the 
weak-coupling expansion! 87 ' 88 The length scale at which 
a weak-coupling expansion breaks down is given by e 7 " 
irrespective of s being integer or half-integeri 87 i 88 For in- 
teger spin chains this implies a gap proportional to the 
inverse cut-off length e _7rs . For half-integer spin chains 
the topological term leads to a different physical behavior 
which resembles that of the s = 1/2 chain, i.e., a gapless 
critical behavior. While this difference is always observed 
in the thermodynamic limit at small fields, it is impor- 
tant to realize that any relevant energy scale such as fields 
or finite-size gaps will lead to a different renormalization 
flow. The physical behavior is then determined by the 
largest energy scale or equivalently the smallest length 
scale. 

In the case of finite chains there are several relevant 
length scales (energy scales), such as the chain length N 
or the correlation length due to finite fields £g in Eq. ([38| . 
The length scale corresponding to the breakdown of the 



weak-coupling expansion is 

N c = e™ (43) 

which corresponds to the correlation length in inte- 
ger spin chains. Two fundamentally different physical 
regimes can be identified (in zero field and temperature): 

(1) N c <C N: This regime corresponds to the most 
studied case of the thermodynamic limit, where the fa- 
mous difference between integer and half-integer spin is 
observed. For finite chains with N > N c (long chains) it 
is possible to clearly see the characteristic features of the 
thermodynamic limit by finite-size scaling (in the form 
of characteristic corrections) . The behavior can be well 
described by continuous quantum field theories; hence 
we call this case the " renormalized continuous quantum 
regime" . 

(2) N <C N c : This regime of short chains was the main 
topic of this paper. The finite-size effects dominate and 
the physical behavior is sensitive to the boundary condi- 
tion and the geometry of the finite cluster, which leads 
to the even-odd effect. It is fundamentally impossible to 
connect the unique behavior in this regime analytically 
to the thermodynamic limit by finite-size scaling. Since 
many of the features are correctly reproduced by the cor- 
responding classical model but quantum effects are still 
important (see below), we call this case the "bounded 
quantum-classical regime" . 

In finite magnetic fields the correlation length £g 
comes also into play. The above two regimes are present 
at low fields, N c <C £,b or N <C £b- The field gives rise 
however to a further regime where £b <C N, N c , which 
we call the " ferromagnetic regime" . It is dominated by 
a relatively large magnetic field or large magnetization 
S ^> 2s, Ns/N c , and both short and long chains enter 
it under these conditions, obliterating the distinction be- 
tween the two low-field regimes. Our numerical results 
show that the correlations are dominated by the trend 
to align all spins with the total spin. This behavior is 
continuously connected to the ferromagnetic region, and 
can be best described by a hydrodynamic theory or by 
spin waves, which give analogous results i 80 ' 89 The behav- 
ior shows no fundamental difference between even and 
odd N nor between integer and half-integer s. For com- 
pleteness it is mentioned that in addition there is also a 
finite temperature regime with smallest length scale J/T, 
which is however not considered in this paper. 

The above considerations demonstrate that the 
crossover from short to long is characterized by the length 
N c , which depends on s and thus describes a boundary 
line as sketched in Fig. [TH] This plot was not completely 
traced out by the numerical results (for obvious limi- 
tations in computational power), but the present work 
provides strong pieces of evidence for its validity. For 
instance, since iV c (l/2) ps 5, the region N -C N c of short 
chains is not reachable in s = 1/2 chains, which is per- 
fectly consistent with Fig. [3J where the variation of gAF 
with N is distinctly different for s = 1/2, and the absence 
of the even-odd effect in the s = 1/2 chain as shown in 
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FIG. 16: (Color online) Sketch of the properties of the AFHC 
in the parameter space spanned by spin magnitude s and 
chain length TV in a s-log TV plot. The region of the short 
chains is separated from the region of long chains by the char- 
acteristic length N c (s) = e ns . The short-chain region con- 
nects to the classical limit s — > oo and the long-chain region 
to the thermodynamical limit TV — > oo. The two regimes are 
suggested to exist also for other, and potentially all antiferro- 
magnetic Heisenberg clusters, and based on the characteristic 
properties they are denoted as "bounded quantum-classical" 
and "renormalized continuous quantum" regime. 



Sec. [V] For s = 1 the short-chain region starts to become 
available, but the short-chain behavior may be realized 
only approximately, which is consistent with e.g. Fig. [5] 
For s > 3/2, however, the short-chain region is already 
available for significant chain lengths. This trend with s 
is also consistent with the notion that in the VBS pic- 
ture the half-integer spin part in the VBS wavefunction 
(when present) becomes less and less relevant the larger 
s is, such that the physics is related to the integer spin 
part. Lastly, the analysis of the classical chain model 
showed that for large s the short-chain behavior is in- 
deed present also in chains with large, but finite TV, as 
expected from Eq. (l43l) . In particular, the slope g(S) is 
suppressed for large even TV in the low field limit, though 
only in a small range of fields B < Js/N. Generally, the 
short-chain region becomes more accessible the larger s 
and the smaller TV is, while the long-chain region is ac- 
cessible for relatively small s and large TV. 

The region of small s and large TV or long chains has 
been studied in great detail in the past and the physics 
can be summarized as the renormalized continuous quan- 
tum regime. For short chains, in contrast, it was demon- 
strated in this paper that many of the features, such 
as the spatial fluctuations, are qualitatively explained 
by the corresponding classical model, and that in this 
sense short chains are classical. However, for experimen- 
tally accessible spin magnitudes s quantum fluctuations 
are clearly not negligible, as demonstrated e.g. by the 
slow convergence of the quantum results to the classi- 
cal limit or the significant overlap of the wavefunctions 
with the VBS states. A further example is the superior 
performance of the Hab model over the classical model 



for the odd chains, and it remains surprising that the 
slope g(S) for odd chains follows the prediction of Hab 
well although the local correlations do not. Short chains 
hence show both classical and quantum aspects, depend- 
ing on the feature one is looking at, and in this sense blur 
the distinction between classical and quantum physics. 
Hence we use the notation "bounded quantum-classical 
regime" for this case. Remarkably, in this regime the 
physical behavior becomes largely independent of s and 
TV, i.e., is generic. 

Implications for spin clusters in general. At this 
point our results are of relevance also from a broader, fun- 
damental perspective. A diagram which at first sight is 
very similar to that in Fig. [TBI was proposed a decade ago 
based on studies on antiferromagnetic Heisenberg rings of 
relatively small size TV but with relatively large spin mag- 
nitudes s.— It was in particular demonstrated that the 
L&_E-band picture and the energies and matrix elements 
calculated with Hab become more accurate the larger s 
and the smaller TV isi 31 i 54 That is, the large s-small TV 
region was (erroneously) linked in this and subsequent 
works to the validity of the Hab model, which lead to 
some inconsistencies »ii In the present work, the situation 
is rectified by identifying the (classical) spatial fluctua- 
tions as the characterizing feature in this regime. Both 
the "small" rings and short chains are apparently in the 
bounded quantum-classical regime. However, for rings 
additionally the Hab model is an excellent approxima- 
tion because of the symmetry-induced homogeneous spin 
densities in their L-band states. For chains, in contrast, 
the Hab model is less appropriate. 

The present work thus states more precisely the dis- 
tinctive behavior of small rings anticipated in Ref. [3l| 
and puts it for the chain systems, through a very de- 
tailed numerical and theoretical analysis, on a solid ba- 
sis. The results lend credit to the idea that the two dis- 
tinctive regimes are in fact generic and present not only 
in rings and chains, but also in other and potentially 
many small antiferromagnetic Heisenberg clusters, since 
one could generally expect that spatial fluctuations dom- 
inate over quantum fluctuations in these systems. For 
the calculation of the quantum fluctuations introduced 
by magnetic anisotropy terms powerful theoretical tools 
exist. 90 In contrast, for the effects of the antiferromag- 
netic Heisenberg interactions in relatively small lattices, 
a satisfying quantitative theory which takes into account 
the relevant effects, i.e., treats the (classical) spatial fluc- 
tuations correctly and introduces the quantum effects, 
appears to be missing at the moment. Developing it 
should be an attractive challenge for the future. 

It is finally mentioned that clusters of a dozen 
exchange-coupled spin centers with relatively large spin 
magnitudes s are currently of high experimental rele- 
vance. An abundance of examples is provided by the 
class of molecular nanomagnets, 91 which through syn- 
thetic chemistry has generated hundreds of magnetic 
molecules with different arrangements of the metal cen- 
ters. The number of spin centers ranges from 4 to a cur- 
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rent maximum of 84, and the spin magnitudes are typ- 
ically s = 3/2, 2, and 5/2 for transition metal clusters. 
Furthermore, the emerging held of the artificially engi- 
neered spin clusters can be expected to provide many fur- 
ther attractive experimental systems^— The key find- 
ings in this work should be of fundamental importance 
to a variety of currently studied experimental systems, 
where mesoscopic effects are very important, and very 
promising applications go hand in hand with interesting 
many-body effects. 



For pairwise interactions we introduce the notation 
Tq(kikj\siSj), indicating a many-spin ITO T k (. ..) in 
Eq. (|A1[) with corresponding values ki and kj for the 
single-spin ITOs T ki (si) and T**(s,-) and all h = if 
I =^ i, j. Note that the elementary ITOs are defined as 



Tt( Si ) = Sf, T^( Sl ) 



1 

: 71 



(A3) 



while Tq (sj) is the identity. For a Heisenberg system only 
the ITO representation of Sj • Sj is needed, which is 
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Appendix A: Spin coupled basis 

In view of the SU(2) symmetry of Eq. (fTJ) it is con- 
venient to perform numerical work directly in a basis of 
eigenstates of the total spin operator S. Since this is 
rarely done in the area of quantum spin systems, some 
details shall be given here. The general procedure is given 
by Racah's methods and the irreducible tensor operator 
(ITO) techniques^ 

In this paper the spin coupling scheme was used where 
at first the first two spins are coupled, and then succes- 
sively the next spin to the previous ones: S = ((. . . ((si + 
S2) +S3) + . . .) + sjy). This yields the SU(2) invariant ba- 
sis states |siS2S'i2S35'i23 . . . snS) with intermediate spin 
quantum numbers S12, S123, • ■ •, Sx...n-i- This basis 
is exploited by expressing the Hamiltonian in terms of 
ITOs^ 

The ITO T k (si) of rank k associated to the spin center 
Si has 2k + 1 components Tr(sj) with q = —k,—k + 
1, . . . , k. Coupling ITOs of different rank and different 
spins is generally achieved through 

T q k ({h}Ah}) = [[...[[T k H^)xT k Hs 2 )pxT k ^ S3 )p 
x ...]**-* x T kN (s N )] k (Al) 



where the kj = kx...j have to be populated according to 
the spin coupling scheme and the intermediate spin quan- 
tum numbers. The tensor product of two ITOs thereby 
reads 

lT k *( Si )xT k t( Sj )}« = {hk 3qi q 3 \KQ)T k q ;{s{)T k i{ S] ). 

Qi,Qj 

(A2) 

By repeated application of Eq. (|A2[) the coupling of the 
ITOs in Eq. (|A1[) can be recast into a sum over the prod- 
ucts of single-spin ITOs and Clebsch-Gordon coefficients 
(kikjqiqj\KQ). 



(s 4 .s 3 )--V3T (ll|s, Sj ). 



(A4) 



The parent Hamiltonians for the VBS wavefunctions 
in Eqs. (|40|) and (l4lT) include also higher-order coupling 
terms (s, • Sj) n , which lead to higher order polynomials 
of the ITOs T ki (si) and Tq. (sj) respectively, e.g. (s, • 

s,) 2 = E 91 , 92 (-i) 91+92 ^ 1 (s l )r g 1 2 (s J )rl 9i (s,)rl (?2 (s,). 

These polynomials can be successively reduced by the 
building up principle^ 

T£(s,)T£( S ,) = J2 (kxk 2qi q 2 \k q )[T k -( Si )T^( Si )] k q 



kq 



(A5) 



with 



[T k >( Si )T k *( Si )] k a = (-ir+ k V2kTi{ ki k2 k ) 

q I s s s I 



( S ||T^( Si )|| S >< S ||T^( Si )|| S > 



( S \\T k (s t )\\s 
The reduced matrix elements are given by 



xT k (s t ). (A6) 



(s\\T k ( Sl )\\s)=k! 



(2s + k + 1 



,1 \V2 



2 k (2k)\(2s-k)\ 
The biquadratic term then becomes- 3 



(A7) 



1— 2 2 

(s,. Sj ) 2 = V / 5T °(22|s iSj ) + ^T (ll|s lSj ) + 

(A8) 

For the s — 2 VBS parent Hamiltonian also the (Sj • Sj) 3 
term is needed for which we obtain 

(Si-Sj) 3 = -\/7T °(33|s l s J )-2\/5r (22|s l s J ) 

(3s 4 2 • s 3 2 - s 2 - s 2 + 2) T °(ll|s iSj ) 



s 2 • s 2 



(A9) 



Finally, the VBS parent Hamiltonians in terms of ITOs 
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read: 



10 

y 



k x Fx 

V5T °(22| Sl s 4+1 ) - -^-T °(ll| Si s i+1 ) 



(A10) 

-^T °(33|s iSi+1 ) + l^T °(22|s 4 s l+1 ) 



r u (ll|s A+1 ) + 28 



(All) 



Appendix B: Iterative construction of VBS 
wavefunctions 



VBS wavefunctions can be calculated in the spin cou- 
pled basis using iteration. For a dimer of two spins s the 
VBS wavefunction is known in the spin coupled basis. It 
is degenerate in S, and the maximal S is s. The wave- 
functions are |0) 2 = l s , s ,0) for total spin S — in an 
obvious notation, and go up to [8)2 = \s,s,s) for total 
spin S = s. If it is known what happens when a further 
spin is attached, i.e., if the VBS wavefunction for an open 
chain of 3 spins with spin s is known, then by iteration 
the VBS wavefunction for any s chain of length N can 
be calculated. In each iteration step the corresponding 
basis vectors are extended by one lattice site, i.e. from 

\Si2,.-.,Sx...n-i)n-i to \Si2, ■ ■ ■ , S1...N-1, S)n- The 
method finds the unnormalized coefficients for a VBS 
wavefunction in the spin coupled basis, and the resulting 
wavefunction has therefore to be normalized after com- 
pletion of the iteration. 

For ans = l VBS chain [Figs. HSJc), (d)] the iteration 
reads 



|0)at = \1) N - U 
\1) N = |1)jv-i- 



(Bl) 



\0) N -i 



The basis functions of 1 1) jv an d |0)jv arc orthogonal. The 
number of relevant basis functions for the s = 1 chain 
grows with N like the Fibonacci numbers. For an s = 2 
VBS chain [Figs-HSfg), (h)] one finds 



\0)n 


= \2)n-i 


Wn 


= |2)jv_i 


|2)iv 


= |2)jv_i 



(B2) 



1)n-i, 
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49" 



|1>jv-i + 



;|0) 



N-l, 



and for an s — 3 VBS chain 
|0)jv = |3) w _i, 



(B3) 
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|2>iv 
|3>iv 



|3}at-i 

|3}at-i 

|3) w _i 
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|2)iY-l, 

\2) N -l^ 



IN-l, 



'28, 



10 |2)iv-i + YY|l)iv-i 



\0) N -i 



For the dimerized s — 3/2 valence bond function 
[Figs. [TBTe). (f)] with only one singlet bond between the 
first two spins one obtains 



\0) 2 n = |3/2) 2JV -i, 

|1>2JV = |3/2)2iV-l 



(B4) 



g|l/2) 2 iv-i, 



|2>2JV 
|1/2> 2 JV-1 

|3/2) 2JV _1 



|3/2) 2JV _i + V2|l/2) 2N-1, 



|2) 2 N-2 
|2) 2 JV-2 



2JV-2, 



g|l)2Af-2 + Y g|0) 2 JV-2, 



and for the dimerized s = 3/2 valence bond function with 
two singlet bonds between the first two spins holds 



\0) 2N = |3/2> 2JV -i 



(B5) 



'32, 



|1>2JV = |3/2) 2J v-i-^/-|l/2) 2JV -i, 



|l/2)2W-i = |l)ajv-2) 



'27, 



|3/2) 2 W-1 — W2N-2 — \ — 0) 2 jV-2- 
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Appendix C: Technical details in the derivation of 
the 0(3) NSLM 

In this appendix the 0(3) NSLM is derived in detail. 
With the decomposition (|T5|) the action in (fTSj) is evalu- 
ated. First, expanding the term Js 2 J2i=i ^« ' ^i+i U P 



to order |l/s| 2 : 



n(x l )n(x l+1 ) 

l(Xi)l(x i+ i) 



l 2 ( Xi ) l 2 (x i+1 ) 



2s 2 



2s 2 



+ (-1) 



i+l 



n(xi) 



n(x l+ i) 



(CI) 

i{xi)' 



Differences of the Neel fields can be approximated by 
derivatives which allow to write n(a;j)n(a;,+i) = 1 — 
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Tj[n(xi) - n(x i+ i)] 2 « 1 - ^[e^n^)] 2 . Then up to a 
constant term: 



[9 a; n(x i )] + 



2 , (l(x,) + l(x t+1 )Y 



2.s 2 



(C2) 



,1 2 



J(x i+ i) 
n(xj) n(x i+ i) 



a 2 -(d x ii) 2 



The alternating term requires a careful treatment of 
the boundary conditions. For periodic boundary con- 
ditions and N even the term can be neglected as can be 
seen by writing n(xi) 1( - Xi + 1 ' ) w -2-d x n(x i+ i)l(x i+ i) and 
nixi+i) 1 -^- rj f5 a; n(a; i )l(a; i ) which directly yields 



j s 2 ^(-i) i+i 



,l(x l+ i) l(x 4 
n(Xi) n(x i+ i) 



i=l 

(C3) 

For open boundary conditions additional boundary terms 
remain. Discussion of these terms is left for later and 
periodic boundary conditions are considered now. Intro- 
ducing L = Na and taking the continuum limit: 



Js 2 fli ■ f2i+i — ► J dx 

i=i J o 



as 2 ( dn 

Ox 



(C4) 

Next, the imaginary part in Eq. (fl~5|) is evaluated. The 
Berry phase is antisymmetric under inversion = 
— uj[— fli], hence 



JV 



is w[fl,*] 



dn(xi,T) 



i=i 
rf> 

—i / dr 
Jo 

■\{x u t)} . 



dr 



(C5) 



In Sec. |TV]the first term in Eq. (|C5|) is of "topological" 
significance. The second one in contrast enters in the 
classical equation of motion. 

Up to an additive constant the total action reads 



A 





Jas 2 




I dx 1 dr 


(-) 


'o Jo 


2 


\dx) 



2J j2 

a 



n x 



dn 



1 



A 



top • 



(C6) 



Completing the square the functional integration over 1 
can be performed, giving 1 = ^j(n x <9 T n). Thus, the 
field 1 generates rotations on n. Note that the constraint 
1 • n = is automatically fulfilled. Overall normaliza- 
tion constants are left out. Finally, the 0(3) NSLM is 
generated, with a term as the effective action for the 
Heisenberg chain in the large s limit, see Eq. (|T7| . 

Finally, Eq. (IC3|) is considered for open bound- 
aries. Using n(xi) 1 ^ Xi + 1 ' rj -±d x n{x i+ i)\{x i+ i) and 

n i x i+i)^iT- ~ f d x n(xi)l(xi) the following boundary 
terms are obtained: 



N-l 



js 2 J2(-iy +1 



l(x i+ i) \(xi) 
n(Xi) n(x i+ i)- — 



-as J [<9 x n(0)l(0) + d x n(L)\(L)] N even 
-as J [9 x n(0)l(0) - d x n(L)l(L)} N odd 



(C7) 



These terms have a scaling dimension of one order higher 
than the bulk terms in the action in Eq. (fTT)) . Note that 
they are of the same order as the higher order terms in 
the Euler-Maclaurin sum formula that occur when go- 
ing from the discrete sum to the continuum integral, see 
Eq. (|C4|) . It is assumed that these terms are small. 
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